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CE-QUAL-W2,  Version  3 

by  Scott  A.  Wells  and  Thomas  M.  Cole 


INTRODUCTION:  CE-QUAL-W2  is  a  two-dimensional  water  quality  and  hydrodynamic  code 
supported  by  the  U.S.  Army  Enginer  Research  and  Development  Center  (Cole  and  Buchak  1995). 
The  model  has  been  widely 
applied  to  stratified  surface 
water  systems  such  as  lakes, 
reservoirs,  and  estuaries  and 
computes  water  levels,  hori¬ 
zontal  and  vertical  velocities, 
temperature,  and  21  other 
water  quality  parameters 
(such  as  dissolved  oxygen, 
nutrients,  organic  matter,  al¬ 
gae,  pH,  the  carbonate  cycle, 
bacteria,  and  dissolved  and 
suspended  solids).  A  typical 
model  grid  is  shown  in  Fig¬ 
ure  1  where  the  vertical  axis 
is  aligned  with  gravity. 

This  technical  note  documents  the  development  of  CE-QUAL-W2,  Version  3,  which  incorporates 
sloping  riverine  sections.  Version  3  has  the  capability  of  modeling  entire  river  basins  with  rivers 
and  interconnected  lakes,  reservoirs,  and/or  estuaries.  Four  example  applications  illustrate  the  use 
of  Version  3. 

DEVELOPMENT  RATIONALE:  CE-QUAL-W2  has  been  in  use  for  the  last  two  decades  as  a 
tool  for  water  quality  managers  to  assess  the  impacts  of  management  strategies  on  reservoir,  lake, 
and  estuarine  systems.  A  predominant  feature  of  the  model  is  its  ability  to  compute  the  two- 
dimensional  velocity  field  for  narrow  systems  that  stratify.  In  contrast  with  many  reservoir  models 
that  are  zero-dimensional  with  regards  to  hydrodynamics,  the  ability  to  accurately  simulate  transport 
can  be  as  important  as  the  water  column  kinetics  in  accurately  simulating  water  quality. 

A  limitation  of  Version  2  is  its  inability  to  model  sloping  riverine  waterbodies.  Models  such  as 
WQRSS  (Smith  1978),  HEC-5Q,  and  HSPF  (Donigian  et  al.  [1984),  have  been  developed  for  river 
basin  modeling  but  have  serious  limitations.  One  issue  is  that  the  HEC-5Q  (similar  to  WQRSS) 
and  HSPF  models  incorporate  a  one-dimensional  (1  -D),  longitudinal  river  model  with  a  1-D,  vertical 
reservoir  model  (1-D  for  temperature  and  water  quality  and  zero-dimensional  for  hydrodynamics). 
The  modeler  must  choose  the  location  of  the  transition  from  1  -D  longitudinal  to  1  -D  vertical.  Besides 
the  limitation  of  not  solving  for  the  velocity  field  in  the  stratified  reservoir  system,  any  point  source 
inputs  to  the  reservoir  section  are  spread  over  the  entire  longitudinal  distribution  of  the  reservoir 
layer. 


Figure  1.  Typical  CE-QUAL-W2,  Version  2  model  grid 
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Other  hydraulic  and  water  quality  models  commonly  used  for  unsteady  flow  include  the  1-D 
dynamic  EPA  model  DYNHYD  (Ambrose  et  al.  1988),  used  with  the  multidimensional  water 
quality  model  WASP.  WASP  relies  on  DYNHYD  for  1-D  hydrodynamic  predictions.  If  WASP  is 
used  in  a  multidimensional  schematization,  the  modeler  must  specify  dispersion  coefficients  to 
allow  transport  in  the  vertical  and/or  lateral  directions  or  use  another  hydrodynamic  model  that 
explicitly  includes  these  effects.  In  addition,  the  Corps  model  CE-QUAL-RIV1  (Environmental 
Laboratory  1995),  is  a  1-D  dynamic  flow  and  water  quality  model  used  for  1-D  river  or  stream 
sections.  None  of  these  models  have  the  ability  to  characterize  adequately  the  hydraulics  or  water 
quality  of  deeper  reservoir  systems  or  deep  river  pools  that  stratify. 

In  the  original  development  of  CE-QUAL-W2,  vertical  accelerations  were  considered  negligible 
compared  to  gravity  forces.  This  assumption  led  to  the  hydrostatic  pressure  approximation  for  the 
z-momentum  equation.  In  sloping  channels,  this  assumption  is  not  always  valid  because  vertical 
accelerations  cannot  be  neglected  if  the  z-axis  is  aligned  with  gravity.  In  addition,  the  current 
Version  2  algorithm  does  not  allow  the  upstream  bed  elevation  to  be  above  the  downstream  water 
surface  elevation.  Because  watershed  and  river  basin  modeling  are  becoming  more  important  for 
water  quality  managers,  providing  the  capability  for  CE-QUAL-W2  to  be  used  as  a  complete  tool 
for  river  basin  modeling  is  an  essential  step  in  improving  the  current  state  of  the  art. 

DEVELOPMENT  APPROACH:  There  were  many  approaches  that  could  have  been  imple¬ 
mented  to  incorporate  riverine  branches  within  CE-QUAL-W2.  By  choosing  a  theoretical  basis  for 
the  riverine  branches  that  uses  the  existing  CE-QUAL-W2  two-dimnesional  (2-D)  computational 
scheme  for  hydraulics  and  water  quality,  the  following  benefits  accrued: 

•  Code  updates  in  the  computational  scheme  affected  the  entire  model  rather  than  just  one  of 
the  computational  schemes  for  either  the  riverine  or  the  reservoir  sections,  leading  to  eas¬ 
ier  code  maintenance. 

•  No  changes  were  made  to  the  temperature  or  water  quality  solution  algorithms. 

•  By  using  the  two-dimensional  framework,  the  riverine  branches  had  the  ability  to  predict 
the  velocity  and  water  quality  field  in  two  dimensions  -  this  has  advantages  in  modeling 
the  following  processes:  sediment  deposition  and  scour,  particulate  (algae,  detritus,  sus¬ 
pended  solids)  sedimentation,  and  sediment  flux  processes  as  well  as  making  Manning’s 
friction  factor  stage  invariant  (see  Wells  (1999)). 

•  Since  the  entire  watershed  model  had  the  same  theoretical  basis,  setting  up  branches  and 
interfacing  branches  involved  the  same  process  whether  for  reservoir  or  riverine  sections, 
thus  making  code  maintenance  and  model  setup  easier. 

The  theoretical  approach  was  to  re-derive  the  governing  equations  assuming  that  the  2-D  grid  is 
adjusted  by  the  channel  slope  (Figure  2). 
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THEORETICAL  BASIS:  Details  of  deriving  the  governing  equations  for  CE-QUAL-W2, 
Version  3  for  the  river  basin  model  are  given  in  Wells  (1997).  Table  1  shows  the  governing 
equations  after  lateral  averaging  for  a  channel  slope  of  zero  (original  model  formulation)  and  for 
an  arbitrary  channel  slope. 


Table  1 

Comparison  of  Governing  Equations  for  CE-QUAL-W2  with  and  without 
Channel  Slope 

Changes  in  Governing  Equations 


Equation 

CE-QUAL-W2  Version  2 

CE-QUAL-W2  Version  3 
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Note:  U,W:  horizontal  and  vertical  velocity,  B:  channel  width,  R.  pressure,  g:  acceleration  due  to 
gravity,  -t  x,x  *:  lateral  average  shear  stress  in  xand  z,  p  :  density,  r\ :  water  surface,  a  :  channel  angle, 
Ux.  x-component  of  velocity  from  side  branch,  q.  lateral  inflow  per  unit  length 
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Numerous  algorithmic  changes  were  made  in  the  CE-QUAL-W2  model.  In  addition  to  the  general 
channel  sloping  feature,  these  changes  include: 

a.  Ability  to  choose  the  following: 

•  Turbulence  closure  models  for  each  waterbody  using  eddy-viscosity  mixing  length 
models. 

•  Varying  vertical  grids  between  waterbodies. 

•  Chezy  or  Manning’s  friction  factor. 

•  Reaeration  formulae  based  on  the  riverine  or  reservoir/lake  or  estuary  character  of  the 
waterbody  or  user-defined  formulations. 

•  Evaporation  models  based  on  theory  or  user-defined  formulations 

b.  Ability  to  linearly  link  one  branch  with  another  or  specify  an  internal  dam  or  internal 
hydraulic  structure(s)  (spillways,  gates,  weirs,  and  pipes)  within  or  between  water  bodies 
(the  pipes  algorithm  is  an  unsteady  1-D  hydrodynamic  sub-model  to  the  core  W2  code 
from  Berger  and  Wells  (1999). 

c.  Effect  of  hydraulic  structures  on  gas  transfer  and  total  dissolved  gas  transport. 

d.  Conservation  of  longitudinal  momentum  at  intersections  between  main  branches  and  side 
branches. 

e.  Effect  of  lateral  inflows  from  tributaries  or  the  lateral  component  of  inflows  from  branch 
intersections  on  the  vertical  eddy  viscosity. 

f  QUICKEST/ULTIMATE  numerical  transport  scheme. 

g.  Implicit  eddy  viscosity  formulation  that  removes  its  timestep  requirement. 

h.  Multiple  user-defined  algal  groups. 

i.  Multiple  user-defined  organic  matter  groups. 

j.  Sediment  diagenesis  model. 

TEST  APPLICATIONS 

Richard  B.  Russell  (RBR)/J.  Strom  Thurmond  (JST).  RBR  and  JST  (previously  known  as 
Clarks  Hill  Lake  as  shown  in  Figure  3)  are  Corps  reservoirs  located  on  the  Savannah  River  between 
Georgia  and  South  Carolina.  JST  is  located  immediately  below  RBR.  The  study  involved  investi¬ 
gating  the  effects  of  proposed  pump-storage  operations  on  temperature  and  water  quality  for  the 
two  reservoirs.  This  required  a  rewrite  of  the  code  in  order  to  dynamically  link  the  two  systems. 
The  linkage  algorithm  was  sufficiently  generalized  to  allow  dynamic  linkage  to  any  number  of 
“waterbodies.”  During  algorithm  development,  it  was  recognized  that  a  natural  extension  would 
be  to  allow  modeling  of  the  river  reaches  between  reservoirs,  thus  providing  a  state-of-the-art 
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Figure  3.  Richard  B.  Russell  and  J.  Strom  Thurmond  (Clarks  Hill)  reservoirs 


waterbasin  management  tool  that  would  handle  reservoir  hydrodynamics,  temperatures,  and  water 
quality  in  a  realistic  manner. 

Prior  to  1 996,  both  reservoirs  were  operated  for  peaking  hydropower.  During  the  course  of  the 
study,  extensive  pump-storage  operations  were  conducted  in  1996  and  data  were  collected  that 
allowed  evaluation  of  the  model’s  ability  to  simulate  the  effects  of  pump  storage  on  temperature 
and  dissolved  oxygen  (DO)  in  the  two  reservoirs. 
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All  hydrodynamic/temperature  calibration  parameters,  including  longitudinal  eddy  viscosity/ 
diffusivity,  bottom  friction,  sediment  heat  exchange,  and  short  wave  solar  radiation  absorption/ 
extinction,  were  set  to  their  default  values  with  the  exception  of  the  wind-sheltering  coefficient, 
which  was  set  to  0.9  for  RBR  and  1.0  for  JST.  The  wind-sheltering  coefficient  is  used  to  multiply 
observed  winds  in  order  to  decrease  the  effective  wind  reaching  the  reservoir  surface. 

Temperature  simulations.  The  model  was  calibrated  to  data  collected  during  1988, 1994,and  1996. 
The  year  1988  represents  a  low-flow  year,  1994  represents  a  high-flow  year,  and  1996  represents 
an  average  flow  year  with  extensive  pump-storage  operations. 

Figures  4-6  show  the  results  of  the  model  simulations  for  the  calibration  years  for  RBR,  while  Fig¬ 
ures  7-9  show  the  results  for  JST.  The  x’s  represent  observed  data  and  the  lines  represent  model 
predictions.  The  x’s  are  scaled  to  represent  ±  0.5  °C.  Additionally,  the  absolute  mean  error  (AME) 
and  the  root  mean  square  error  (RMS)  are  presented  on  the  plots  as  a  further  aid  in  evaluating  model 
predictions. 

The  results  show  that  the  model  is  capable  of  reproducing  the  observed  temperatures  in  both  RBR 
and  JST  under  varying  flow  years  and  reservoir  operations.  Additionally,  the  model  captured  the 
approximate  4  °C  increase  in  hypolimnetic  temperatures  for  1996  compared  to  1988  and  1994  due 
to  pump-storage  operations. 
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Figure  5.  1 994  RBR  computed  versus  observed  temperature 
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Figure  7.  1 988  JST  computed  versus  observed  temperature 


Figure  8.  1 994  JST  computed  versus  observed  temperature 
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Figure  9.  1 996  JST  computed  versus  observed  temperature 


Dissolved  oxygen  simulations.  An  additional  complication  to  this  study  was  the  presence  of  an 
oxygen  diffuser  system  that  had  to  be  incorporated  into  the  model  in  order  to  reproduce  observed 
DO  in  RBR.  The  oxygen  diffuser  system  consists  of  one  system  located  in  the  forebay  and  another 
located  approximately  1.6  km  upstream  from  the  dam.  During  1988,  most  oxygen  injection 
occurred  in  the  forebay  except  during  late  summer  when  the  upstream  system  was  also  used. 
Oxygen  injection  during  1994  utilized  both  of  the  systems,  while  in  1996  the  upstream  system  was 
the  only  one  utilized. 

Figures  10-13  show  the  results  of  1988,  1994,  and  1996  DO  simulations  for  RBR  for  the  station 
located  closest  to  the  dam.  The  simulations  included  the  full  suite  of  algal/nutrient/DO  interactions 
available  in  the  model.  Sediment  oxygen  demand  (SOD)  was  modeled  using  a  first-order  sediment 
compartment  that  included  the  effects  of  organic  matter  delivery  to  and  subsequent  decay  in  the 
sediments  for  both  allochthonous  and  autochthonous  sources  of  organic  matter.  Additionally,  a 
zero-order  sediment  oxygen  demand  was  included  to  provide  for  background  SOD. 

Observed  DO  in  RBR  exhibits  very  complex  patterns  in  response  to  the  oxygen  injection  system, 
algal  production,  peaking  hydropower  operations,  and,  in  1 996,  pump-storage  operations.  For  1988 
and  1 994,  the  DO  regime  exhibits  a  pronounced  metalimnetic  minimum  due  to  the  oxygen  injection 
system  delivering  oxygen  to  the  hypolimnion.  In  contrast,  the  metalimnetic  DO  minimum  is  greatly 
reduced  in  1996  due  to  the  effects  of  pump  storage  that  caused  greater  hypolimnetic  mixing 
(Figure  13).  Although  not  normally  thought  of  when  calibrating  for  hydrodynamics,  DO  in  RBR 
provides  a  much  better  indication  of  how  well  the  model  is  reproducing  reservoir  hydrodynamics 
because  of  the  very  dynamic  interactions  of  peaking  hydropower,  pump  storage,  and  oxygen 
injection  operations. 
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Figure  11.  1988  RBR  computed  versus  observed  DO,  June-October 
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Figure  1 3.  1 996  RBR  computed  versus  observed  DO 
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Figure  14  illustrates  the  model’s  ability  to  capture  the  longitudinal  gradients  in  DO  along  the 
mainstem  as  well  as  the  differences  in  DO  between  side  branches  and  the  mainstem  during  July  of 
1996.  CE-QUAL-W2  reproduces  the  anoxic  hypolimnion  for  the  two  stations  located  off  the 
mainstem  and  reproduces  the  differences  between  the  two  stations.  The  model  also  reproduces  the 
effects  of  additional  hypolimnetic  oxygen  delivery  upstream  of  the  oxygen  diffuser  system  due  to 
pump-storage  operations  advecting  higher  DO  upstream. 


Figure  1 4.  July  1 996  RBR  computed  versus  observed  DO.  The  first  two  plots  are  stations  located  on  side 
branches  off  the  mainstem.  The  last  five  plots  are  stations  located  at  the  upstream  portion  of 
the  mainstem  down  to  the  dam 


Figures  15-17  show  the  results  of  the  1988, 1994,  and  1996  DO  simulations  for  JST.  The  DO  regime 
in  JST  is  quite  different  from  that  in  RBR,  although  JST  also  exhibits  a  metalimnetic  DO  minimum. 
The  model  captures  much  of  the  DO  dynamics  in  JST.  However,  the  metalimnetic  DO  minimum 
is  not  captured  as  accurately  as  in  RBR. 
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Figure  1 5.  1 988  JST  computed  versus  observed  DO 
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Figure  16.  1994  JST  computed  versus  observed  DO 
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Figure  1 7.  1 996  JST  computed  versus  observed  DO 


Lower  Snake  River.  The  domain  of  the  Lower  Snake  River  from  C.  J.  Strike  Reservoir  (RM 
487)  to  the  headwaters  of  Brownlee  Reservoir  (RM  335)  is  244  km  (152  miles)  in  length.  The  river 
was  broken  into  five  branches  of  varying  slope  from  0.001  to  0.0008  (Figure  18).  The  model 
consisted  of  312  longitudinal  segments  between  805  and  835  m  in  length,  13  tributary  and  point 
sources,  1  distributed  load,  and  90  agricultural  return  flows. 

Hydraulics  were  calibrated  using  water  surface  elevation  data  at  specific  flow  rates.  Gauging  station 
data  were  available  at  several  locations  throughout  the  domain.  Figure  19  shows  the  water  level 
calibration  for  a  flow  of  5,600  cfs.  Mean  water  level  error  and  root  mean  square  water  level  error 
for  flow  rates  between  5,600  cfs  and  13,000  cfs  were  well  below  0.5  ft  for  a  river  that  experiences 
a  300-ft  drop  over  its  length.  The  calibrated  Chezy  values  varied  from  segment  to  segment  between 
20  and  80  and  were  flow  and  stage  invariant. 

The  primary  goal  of  this  modeling  study  was  to  determine  the  loading  of  organic  matter  and  nutrients 
to  Brownlee  Reservoir.  Model  predictions  of  temperature,  algae,  nutrients,  and  organic  matter 
compared  well  with  field  data  at  six  locations  along  the  river.  Results  from  the  Porter’s  Island  station 
are  given  in  Figures  20-23. 

Bull  Run  River  System.  The  Bull  Run  watershed  has  been  the  primary  drinking  water  supply 
since  1 895  for  the  metropolitan  area  of  Portland,  OR,  USA.  The  watershed  is  composed  of  two 
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Figure  1 8.  Lower  Snake  River  model  segmentation 


Figure  1 9.  Comparison  of  model  predictions  and  data  of  water  surface  along  a  portion  of  the  Lower 
Snake  River  at  a  flow  of  5600  cfs 
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Figure  20.  Computed  versus  observed  temperatures 


Figure  21 .  Computed  versus  observed  chlorophyll  a 
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Figure  23.  Computed  versus  observed  nitrate/nitrite 
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man-made  reservoirs  (Reservoirs  1  and  2),  and  a  potential  third  reservoir.  Because  of  compliance 
requirements  for  endangered  species  survival,  the  reservoirs  and  river  segments  in  the  watershed 
were  modeled  with  Version  3  in  order  to  meet  temperature  standards  for  fish.  Figure  24  shows  the 
profile  of  the  model  system,  including  two  river  and  three  reservoir  sections.  River  channel  slope 
was,  on  average,  greater  than  2  percent  in  the  Upper  Bull  Run  River. 
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Figure  24.  Bull  Run  River-Reservoir  system  profile 

A  dye  study  performed  in  the  Lower  Bull  Run  River  during  June  1 999  verified  the  river  model  travel 
times  and  dispersive  characteristics.  Model  results  and  observations  are  compared  in  Figure  25 
using  the  QUICKEST/ULTIMATE  numerical  scheme. 

Model  predictions  of  temperature  profiles  in  the  two  reservoirs  during  a  two-year  continuous 
simulation  period  were  within  0.5  °C  Absolute  Mean  Error  and  0.6  °C  Root  Mean  Square  error  for 
over  40  profile  comparisons  in  each  reservoir.  A  typical  series  of  model-data  predictions  for 
Reservoirs  1  and  2  are  shown  in  Figures  26-29.  Note  how  the  model  captures  the  double  thermocline 
in  late  summer  for  Reservoir  1  and  the  differences  between  the  thermal  regimes  in  the  two  reservoirs. 
All  hydrodynamic/temperature  coefficients  were  set  to  their  default  values  except  for  wind- 
sheltering. 
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Figure  25.  June  1 999  Lower  Bull  Run  River  computed  versus  observed  dye  concentrations 
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Figure  26.  1 997  Bull  Run  Reservoir  1  computed  versus  observed  temperatures 
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Figure  27.  1 998  Bull  Run  Reservoir  1  computed  versus  observed  temperatures 


Figure  28.  1 997  Bull  Run  Reservoir  2  computed  versus  observed  temperatures 
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Figure  29.  1 998  Bull  Run  Reservoir  2  computed  versus  observed  temperatures 

Columbia  Slough  System.  The  Columbia  Slough  (Figure  30)  is  an  extensive  system  of 
interconnected  wetlands,  channels,  and  lakes  located  in  the  Portland,  Oregon,  USA  metropolitan 
area  and  lying  in  the  floodplain  of  the  Columbia  River.  It  is  approximately  30  km  in  length  and 
includes  a  freshwater  estuary  portion  and  a  series  of  isolated  lakes  and  channels  that  receive 
stormwater  and  groundwater  inflows.  The  model  was  developed  to  evaluate  the  effect  of  combined 
sewer  overflows,  stormwater,  and  groundwater  inflows  on  water  quality  in  the  Columbia  Slough 
system.  Model  development  is  summarized  in  Berger  and  Wells  (1999). 

The  Lower  Columbia  Slough  (Figure  3 1 )  is  connected  to  the  Willamette  River,  where  it  experiences 
a  water  surface  fluctuation  between  1  and  3  ft.  Inflows  to  the  Lower  Columbia  Slough  include 
combined-sewer-overflows  (CSOs),  storm  water,  water  from  Smith  and  Bybee  Lakes,  leachate  from 
the  St.  John’s  Landfill,  and  inflows  (both  pumped  and  gravity  inflows)  from  the  Upper  Columbia 
Slough  at  MCDD1. 

The  Upper  Columbia  Slough,  shown  in  Figure  32,  was  historically  maintained  to  provide  irrigation 
water  to  agricultural  and  commercial  users  in  the  summer  months.  The  Upper  Slough  is  connected 
by  pipes  and  an  overflow  weir  during  the  fall,  winter,  and  spring  to  Fairview  Lake.  During  the 
summer,  Fairview  Lake  is  connected  to  the  Upper  Slough  only  by  flow  over  and  leakage  through 
the  weir.  Water  is  also  pumped  from  the  Upper  Slough  to  the  Lower  Columbia  Slough  at  a  pump 
station  at  MCDD1  (Multnomah  County  Drainage  District  No.  1)  and  from  the  Upper  Slough  to  the 
Columbia  River  at  MCDD4  (Multnomah  County  Drainage  District  No.  4).  At  MCDD1 ,  pipes  also 
allow  gravity  flow  from  the  Upper  Slough  to  the  Lower  Slough.  Other  inflows  to  the  Upper 
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Figure  30.  Columbia  Slough 


Figure  31 .  Lower  Columbia  Slough 
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Figure  32.  Lower  Slough  measured  center-line  vertical  velocity  profiles  compared  to  laterally  averaged 
model  predictions  during  high-water  conditions. 


Columbia  Slough  include  groundwater  and  storm  water  from  the  Portland  International  Airport  and 
other  industrial,  commercial,  and  residential  neighborhoods  in  the  area. 

The  model’s  ability  to  reproduce  velocities  in  the  tidally  dominated  Lower  Slough  is  shown  in  Fig¬ 
ure  33  during  high-water  conditions.  CE-QUAL-W2  velocity  predictions  are  laterally  averaged, 
whereas  velocity  measurements  were  taken  at  the  channel  center. 

The  model’s  ability  to  reproduce  water  surface  elevations  in  the  Upper  Slough  is  shown  in  Figure  34. 
Note  that  the  model  is  capable  of  reproducing  the  effects  of  water  hammer  on  the  water  surface  due 
to  pumping  operations. 

The  model’s  ability  to  reproduce  flows  in  the  Upper  Slough  is  shown  in  Figure  35. 

SUMMARY:  A  2-D  hydrodynamic  and  water  quality  model  (CE-QUAL-W2,  Version  3)  that 
allows  the  integration  of  river,  reservoir,  lake,  and  estuarine  systems  has  been  developed  for  river 
basin  modeling.  Four  test  cases  demonstrate  the  ability  of  the  model  to  reproduce  river  and  tidal 
hydraulics  and  temperature  and  DO  dynamics  in  stratified  reservoirs. 

Additional  improvements  in  Version  3  are  being  explored  including  the  application  of  a  k— 8 
turbulence  model  rather  than  the  existing  mixing  length  model  for  the  vertical  transfer  of  momentum 
and  an  implicit  solution  for  water  surface  elevations  among  branches.  A  pre/postprocessor  is 
undergoing  development  and  will  be  updated  for  Version  3.  Version  3  is  still  undergoing  algorithm 
development  and  testing.  It  will  be  officially  released  at  the  next  CE-QUAL-W2  workshop.  Details 
on  the  workshop  can  be  found  at  http://www.ce.pdx.edu/~scott/w2 . 
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Figure  33.  Columbia  Slough  computed  versus  observed  water  surface  elevations 
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Figure  35.  Upper  Slough  computed  versus  observed  flows 


POINTS  OF  CONTACT:  This  technical  note  was  written  by  Dr.  Scott  Wells,  Portland  State 
University,  and  Dr.  Thomas  Cole,  U.S.  Army  Engineer  Research  and  Development  Center.  For 
additional  Information,  contact  Dr.  Cole  (601-634-3283,  tcole@wes.army.mil )  or  the  Managers  of 
the  Water  Quality  Research  Program,  Dr.  John  Barko  (601-634-3654,  barkoj@wes.army.mil )  and 
Mr.  Robert  C.  Gunkel  (601-634-3722,  gunkelr@wes.army. mil).  This  technical  note  should  be  cited 
as  follows: 

Wells,  S.  A.,  and  Cole,  T.  M.  (2000).  “CE-QUAL-W2,  Version  3,”  Water  Quality 
Technical  Notes  Collection  (ERDC  WQTN-AM-09),  U.S.  Army  Engineer  Research  and 
Development  Center,  Vicksburg,  MS.  www.wes.army.mil/el/elpubs/wqtncont.html 
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